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ABSTRACT 

We present a Semi-Analytical Line Transfer model, SALT, to study the absorption and re-emission 
line profiles from expanding galactic envelopes. The envelopes are described as a superposition of 
shells with density and velocity varying with the distance from the center. We adopt the Sobolev 
approximation to describe the interaction between the photons escaping from each shell and the 
remaining of the envelope. We include the effect of multiple scatterings within each shell, properly 
accounting for the atomic structure of the scattering ions. We also account for the effect of a finite 
circular aperture on actual observations. For equal geometries and density distributions, our models 
reproduce the main features of the profiles generated with more complicated transfer codes. Also, our 
SALT line profiles nicely reproduce the typical asymmetric resonant absorption line profiles observed 
in starforming/starburst galaxies whereas these absorption profiles cannot be reproduced with thin 
shells moving at a fixed outflow velocity. We show that scattered resonant emission fills in the resonant 
absorption profiles, with a strength that is different for each transition. Observationally, the effect 
of resonant filling depends on both the outflow geometry and the size of the outflow relative to the 
spectroscopic aperture. Neglecting these effects will lead to incorrect values of gas covering fraction 
and column density. When a fluorescent channel is available, the resonant profiles alone cannot be 
used to infer the presence of scattered re-emission. Conversely, the presence of emission lines of 
fluorescent transitions reveals that emission filling cannot be neglected. 

Subject headings: galaxies: ISM ISM: structure 


1. INTRODUCTION 

The mechanical and radiative energy (from massive 
star winds and supernova explosions) injected in the in¬ 
terstellar medium (ISM) of galaxies is expected to drive 
gas outflows around regions of active star-formation. 
These outflows are indeed observed both on the scale 
of individual H II regions as well as on full-galaxy 
scales. Outflows are currently invoked as the princi¬ 
pal mechanism regulating the galactic baryonic cycle 
(i.e., the balance between the gas accretion rate and 
the star-formation rate) in state-of- the -art g alax y forma¬ 
tion models (e.g . iQppenheimer et~ahl 120101 : iDave et all 
120111 ILillv et ahl I2013H . Whether or not the observed 
outflows are actually able to do the job is, however, 
still unclear. In fact, although outflow-regulated models 
are able to broadly reproduce some of the fundamental 
correlations observed in massive galaxies, recent studies 
have pointed out the existence of a fundamental prob¬ 
lem with the evolution of low mass (M* < 10 9 ' 5 M Q ) 
galaxies (iWeinmann et al1l2012r> . This difficulty appears 
to indicate that crucial feedback processes are modeled 
incorrectly, since it is precisely in this low-mass regime 
that feedback-induced outflows are expected to have the 
strongest impact. 

What is lacking at this point are robust observational 
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constraints on the physical properties of galaxy out¬ 
flows and how these depend on the galaxy star-formation 
rates, stellar masses, and so on. In order to characterize 
how effective feedback is in quenching the star-formation 
(by, e.g., heating the gas and/or completely remove it 
from a galaxy dark matter halo) we need to be able to 
probe the kinematics of the outflows, their density struc¬ 
ture, and their extent. Absorption line studies against 
the strong UV continuum produced in the star-forming 
regions can probe the neutral and ionized components of 
the outflows, with numerous resonant transitions in low 
and high ionization metals ( such as, Si II, C II , Mg II, 
C III , Fe II, andC IV; e.g.. iRupke et all 120051: iMartin 
2005 : ISato et al.1 l2009t iWeiner et ah 20091 Rubin et al 


201C: ISteidel et ah 1996; Shaplev et al JDOOdfT 


Recent studies show that outflowing gas moving with 
velocities up to several hund reds of km s" 1 is common 
in sta r-forming galaxies (e.g.- IPettini et alll2002t[Martini 
I2005T) . Generall y, these studies use standa rd absorption¬ 
line analysis (e.g lSavage fc Sembachll 19911) . where the ve¬ 
locity of the outflowing material is determined by the 
amount of blueshift observed in the resonant absorption 
lines (typically UV lines, but absorption of Na- D lines 
has been used). Recent observations, however, have re¬ 
vealed the presence of numerous resonant and fluorescent 
emission lines associated with the blueshifted resonant 


absorptions ( 

e.g.. IWeiner et al.l 

20091: iFrance et al.l 

2010; 

Rubin et al.l 

2011 

Jones et al 

201 § Erb et alJ 

2012; 

Martin et al.l 

2012 

20131: IKornei et al. 201311. Although 


originally interpr eted as the res ult of photoionization by 
weak AGNs (e.g.. IWeiner et al1l2009fl . these lines are now 
believed to be the result of scattered reso nant photons 
in an expanding envelope around galaxies dRubin et al.l 
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Fig. 1.— Side view of the outflow model. We consider a spherical galaxy (with radius Rsf) surrounded by an expanding envelope with 
a velocity increasing radially with the distance from the center. For r > Rw the gas velocity is Uoo- The shaded yellow area shows the size 
of the aperture of radius R a per- The dotted vertical line indicates the edge-on view of the plane with constant observed velocity (v 0 b s )- 


[20ntlErb et~ahll20lllMartin et al.ll2013lf . 

The scattered re-emission into the line of sight affects 
the velocity measurements based on pure absorption lines 
analysis, as well as mi mic a partial covering f raction of 
the outflowing gas le.g.. iProchaska et alJl2011[ ). In these 
cases, it is crucial to be able to consistently model both 
the resonant absorption as well as the associated res¬ 
onant and (if present) fluorescent emission originating 
from the same io nic species. Steps in this direction have 
been taken by iRubin et al.1 (120Ilf) , who pioneered the 
study of outflows in emission using resonantly-scattered 
Mgll, and Fe lines. Prochaska et al. (2011) used Monte 
Carlo radiative transfer techniques to study the nature of 
resonant absorption and emission for winds, accounting 
for the effects of resonant scattering and fluorescence. 

In this paper we go a step forward and develop a Semi- 
Analytical Line Transfer (hereafter referred to as SALT) 
model to interpret the absorption/scattered emission line 
profiles resulting from extended galactic outflows. We 
assume that the Sobolev approximation holds, and we 
account for multiple scattering within the outflow with 
a simple statistical approach. As an example of an ap¬ 
plication of the model, we apply SALT to multiple tran¬ 
sitions in the Si + ion observed in the stack spectrum of 


z ~ 0.3 Lya emitters. We show how our model is able to 
consistently reproduce the profiles of both the absorption 
and resonant and fluorescent emission lines and thus con¬ 
strain the outflow velocity- and density-fields. Although 
rather simplified, the SALT model can be used to gain 
a more physical understanding of the outflowing gas, in 
both local and high redshift galaxies. We focus here on 
a few specific lines of Si + ion, however, the models are 
readily applicable to any other ions. 

The paper is organized as follows. In Section 2 we 
present the derivation of the semi analytical line profile 
for an outflowing spherical shell. The model is compared 
with the observed stacked spectrum in Section 3, and the 
results are discussed in Section 4. We offer our conclu¬ 
sions in Section 5. 

2. P-CYGNI PROFILE FROM A SPHERICAL EXPANDING 
ENVELOPE 

In spherical outflows, such as those produced in winds 
from early type stars, resonant lines are chara cterized 
by the well understood P-C ygni profile (e.g., iCastorl 
119701 : iCastor fe Lamerslll979D . The profile shows both 
an emission and an absorption component. The absorp¬ 
tion is created in the material between the source and 
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the observer, while the emission is produced by scatter¬ 
ing photons into the line of sight. For spherical outflows, 
the velocity profile of the absorption component will be 
blueshifted relative to the systemic velocity of the source 
and will depend on the density and velocity fields of the 
absorbing material. Because of the spherically symmet¬ 
ric geometry, the emission component will be centered 
at the systemic velocity of the source and thus will con¬ 
tribute to fill in the absorption at negative velocities. 
This effect, well known in the context of stellar winds, 
has generally been neglected in the context of absorp¬ 
tion line studies of galaxies. Recently, however, a few 
studies have emphasized the importance of properly ac¬ 
counting for scattered re-emission in studies of galaxy 
absorption line spectra. These studies have also high¬ 
lighted how the emission line features can be used as 
powerful diagnostics of the geometry and physical con¬ 
ditions of galaxy o utflows both in th e local and high- 
redsh i ft Universe (iRubin et al.l 120111: iProchaska et al.l 
1201 It lErb et all l2012t iMartin et alJ 12013ft .The goal of 
this paper is to present a semi-analytical line transfer 
(SALT) model that can be used to consistently interpret 
the absorption and emission line profiles of both reso¬ 
nant and fluorescent transmissions observed in galactic 
spectra. The line profiles are modeled for outflow ge- 
o metries similar to those introduced by the recent works 
oflRubin et al.1 (12011 1: IProchaska et al.l (1201 If ) : lErb et al.1 
(120121k iMartin et, ail (|2013| h 
In this Section we first describe the basic assumptions 
to deri ve the lin e pr ofiles in sing le s cattering approxi¬ 
mation |ScuderT^OD (following Il992f) . We then modify 
the wind model to include a more realistic description 
of the scattering process by relaxing the single scatter¬ 
ing approximation. We also include the possibility of 
re-emission in the fluorescent channel, different envelope 
geometries and effect of a finite aperture. 

2 .1. The basic model 


introduce the normalized radial coordinate g = t/Rsf 

(g = \/C 2 + s 2 /Rsf)- A given point P in the envelope is 
identified by the pair of coordinates p, 9 1 where 9 is the 
angle between the direction of r and the line of sight to 
the observer. 


Face on view 


6 ... 1.0 
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Fig. 2.— Iso-velocity contours of observed velocity for a shell 
moving at radial velocity v. The observed velocity ranges from 
u obs = v at the center of the shell, to u obs = 0 at the edge. 

We consider a velocity field where the velocity (v) in¬ 
creases with r as a power law of exponent 7 : 

u = for r - R w = / ; (1) 

v = Vqo for r > R w - (2) 


To build our SALT model we start with the simplest 
description for the gas/star configuration. We approxi¬ 
mate a galaxy as a spherical source of UV radiation with 
radius Rsf (i.e., where the bulk of star-formation oc¬ 
curs), surrounded by an expanding envelope of gas, ex¬ 
tending to R\v- In the following we will refer to the ex¬ 
panding envelope as the “galactic wind”. Prochaska et al. 
(2011) use a Monte Carlo radiative transfer technique to 
derive the absorption line profiles resulting from similar 
galactic winds. In what follows, we simplify their calcu¬ 
lations by computing semi-analytical expressions for the 
line profiles (including the effects of fluorescent emission 
and spectroscopic aperture size). Martin et al. (2013) 
use a similar outflow model to study the extended scat¬ 
tered emission from galactic outflows at z ~ 1, consid¬ 
ering the effects of a clumpy expanding medium on the 
derived mass outflow rate. Here, we build upon these 
works, and we develop a semi-analytical algorithm to 
simply but accurately calculate the expected line pro¬ 
files originating in extended outflows. Our results can 
be used to easily model galactic extended outflows com¬ 
monly observed in local and high redshift galaxies. 

In Figure [T] we show the geometry of the wind and the 
definition of the coordinate system. The coordinates £ 
and s are given in units of Rsf , so that in Figure [T| the 
dashed line tangential to the galaxy has £ = 1. We also 


where Vq is the wind velocity at the surface of the star¬ 
forming region (i.e., at Rsf), and Voo is the terminal 
velocity of the wind at Rw- 
When the velocity gradient in the expanding envelope 
is large, photons will interact with the outflowing mate¬ 
rial only where the absorbing ions are exactly “at reso¬ 
nance” due to their Doppler shift (this is a condition als o 
known as “Sobolev’ approximation”, e.g., iGrini nl l200l . 
In this case, the radiative transport of the line photons 
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Fig. 3.— Energy levels of Si+. 
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can be reduced to a local problem, and the optical depth 
for absorption (t) can be evaluated at the interaction 
surface , which is defined in terms of the velocity as: 


Au 

’ 


( 3 ) 


where uq is the resonance frequency of the line. The 
numerical results of Prochaska et al (2011) show that for 
the physical conditions of typical galaxy outflows, the 
Sobolev approximation is a justified assumption. 

The wind optical depth at the interaction surface, can 
be written as a function of wavelength, wind parame - 
ters, and atomic constants, as follows ('e.g.. lCastodll970f) : 


r(r) 


7re 


fiuM u ni{r) 


n u gn r/v 
nig u \ 1 + erg 2 ’ 


( 4 ) 


where f u i and A i u are the oscillator strength and wave¬ 
length, respectively, for the ul transition, g = cos( 0 ), 

and er = — 1 (see Table Q] for Si + atomic data). 

The expression for the optical depth can be simplified 
by assuming that 1 ) it does not depend on the angle 
9 between a radius and the line of sight, 2) the above 
velocity law holds, 3) stimulated emission is negligible 


(i.e., 


nufll 

nig u 


1), and 4) the mass outflow rate is 


constant, so that rq(r) oc (it 2 ) . For 7 = 1, and with 

these assumptions, we can write: 


and 


( \ ^ t \ ( Rsf \ r 

r(r) = - JiuMuTio - 

me \ r J v 

(RsfV 
= T0 {—) 


( 5 ) 

( 6 ) 


TO flu^lu^l 0 ? (0 

me vo 

where no is the gas density at Rsf (for 7 = 1 , ni(r) = 

»»(7fe)' 3> ' 

Now consider a thin shell located at a distance r = 
Rsf(^) 1 7 , moving with an intrinsic radial velocity v. 
The velocity measured by the observer, i.e. the compo¬ 
nent of the radial velocity along the line of sight to the 
observer (r> 0 bs = v cos 9) will depend on the position on 
the shell and in particular on the projected distance to 
the center. This is shown in Figured where we plot con¬ 
tours of constant observed (i.e., projected) velocity from 
a shell moving outward with radial velocity v. For the 
sake of clarity, we show only the half of the shell moving 
toward the observer. The observed velocities range from 
—v at the projected center of the envelope (where the 
gas is moving directly toward us) to 0 at the projected 
distance r = r v (where the shell is moving on the plane 
of the sky). Obviously, only the portion of the shell in 
front of the continuum disk (hatched area in Figures 0 
will produce a net absorption in the spectrum (blueward 
of the line center) by scattering photons out of the line 
of sight. 


Globally the shell will absorb a fraction E(v) = [1 — 
exp (—t(u))] of the energy that, in terms of observed ve¬ 
locities (fobs = v cos 9) , will be redistributed evenly over 
the velocity interval (u m in, v). Here i> m i n is the projec¬ 
tion of the shell velocity along the line of sight tangential 
to the galaxy (i.e., at £ = 1). Following Scuderi et al. 
(1992), we can compute v m i n as: 

v min =v cos9 = v^y^-. ( 8 ) 

e{v) 

Or, setting y = v/vq, as: 

2 /min = y (7-1)/7 (y 2h - i) 1/2 . (9) 

Only shells with intrinsic radial velocities in the range 
from i> 0 b s and v± = nobs/ cos 9 (for £ = 1) can contribute 
to the absorption at n Q b s - Setting x = n 0 bs/no, yi = 
Vi/vq can be computed by solving the equation: 

yl{l~yf 2,1 )=x 2 . (10) 

Thus, we can write the absorption component of the 
profile, in units of the stellar continuum as: 

ryr i _ e ~ T (y) 

-Gbs,blue(*n) — / d y. (H) 

J max{ x,l) 9 2/min 

Surfaces of constant observed velocity can be described 
by the equation: 


^obs — ^0 



( 12 ) 


For the particular case of 7 = 1, this equation describes 
parallel planes at distance r = RsF'^f- cos 9 from the 
center of the emitting region (see Figure [1]). 


2.1.1. Single scattering approximation 

Resonant photons absorbed in the envelope can be de¬ 
tected when re-emitted toward the observer. Assuming 
that a re-emitted photon escapes a given shell without 
further interactions, we can compute the emission com¬ 
ponent of the line profile as follows. If the photons are re¬ 
emitted isotropically, then they will uniformly cover the 
range of projected velocities between ±v. To describe the 
profile of this emission component, we divide the range 
of observed velocities into blueward and redward of the 
systemic velocity. The blue side of the emission profile 
originates in the half of the envelope approaching the 
observer (i.e., s > 0). For a given observed velocity, the 
emission will come from all shells with v > u 0 bsi an< 4 we 
can write: 


f v °° 1 — e~ T ^ 

4em,blue(^) = / ^ dy. (13) 

Jmax(x, 1) ^V 

The red side of the profile is produced in the unoc¬ 
culted, receding portion of the envelope. Because of the 
occultation, however, only shells with velocities larger 
than v m i n (see Eq. [9]) will contribute to a given observed 
receding velocity: 


7em,red (-C ) 



1 - e~ T ( y ) 

2 y 


d y. 


(14) 



















5 


Finally, the resulting P-Cyg profile for the ideal spher¬ 
ical outflows can be computed as: 

I(%) “ 1 F a bs,blue T ^em,blue T d em? red- (15) 


outflow, where the density is highest. We account for 
multiple scatterings within a single shell, as follows. 

We define the photon’s esca pe prob ability from a shell 
of optical depth t(v) as fe.g.. lMeithislll972 d: 


2.1.2. Fluorescent emission in single scattering 
approximation 

Depending on the energy levels of the particular ionic 
species, the absorption of a resonant photon can result 
in the production of a fluorescent photon. This occurs 
when the electron decays into an excited ground leve@. 
As an example, figure [3] shows the energy level diagram 
of the A 1190.42 and 1193.28A Si + doublet. Resonant 
and fluorescent transitions are marked with dashed and 
dotted lines, respectively. We account for the fluorescent 
channel in the modeling of the line profile as follows. 

For a bound electron, the probability of decaying into 
the lower level l is proportional to p u i = A u i/J2i A U i, 
where A ui is the spontaneous decay probability from the 
upper level u to the lower level i. Relevant A U i values are 
given in Table [T] In the single scattering approximation, 
the resulting line profile accounting for the emission in 
the fluorescent channel becomes: 


d(x) — 1 A.bs,blued - !?.R (Am,blued - -Am,red)d~P.F (Am,blued - A] 

(16) 

where pr and pf are the probabilities that a photon is 
re-emitted in the resonant and the fluorescent channels, 
respectively. When only one fluorescent channel is avail¬ 
able, as in the cases considered here, pr+pf = 1- The 
left panel in Figure [4] shows how the line profiles of the 
Si II doublet generated in an outflowing envelope change 
when the fluorescent channels are taken into account. 
When a fraction of the photons are re-emitted in the 
fluorescent transition, the filling effect of the resonant 
absorption due to photons scattered in the wind is sub¬ 
stantially reduced for the 1190A transition. It is only 
minimally reduced for the transition at 1193A due to 
contamination from the fluorescent re-emitted photons 


at 1194.5A (see Table [T]). Because of the re-emission in 
the fluorescent channel, the total equivalent width of the 
resonant P-Cygni profile (i.e., including both the absorp¬ 
tion and emission components) is negative (net absorp¬ 
tion). 


2.1.3. Accounting for multiple scatterings 

More realistically, a photon re-emitted with the reso¬ 
nant energy will likely interact with the ions in the shell 
where it was created, resulting in multiple scattering 
events of a single photon within a given shell. Multi¬ 
ple scatterings will not change the shape of the absorp¬ 
tion profile (labs,blue), but will reduce the contribution 
of the re-emission in the resonant line, while enhanc¬ 
ing the re-emission in the fluorescent channel (when this 
is available). The number of scatterings will clearly be 
a function of the ion density at any given point in the 
outflow. Thus, for our assumed density profile, this pro¬ 
cess will be more important in the internal regions of the 


/3 = (l-e t )/t, (17) 

Thus, for a shell with velocity v, a photon has a prob¬ 
ability /3 of escaping the shell, and therefore -because 
of the underlying Sobolev approximation- escaping the 
outflow. Of all photons absorbed at resonance by the 
moving shell, a fraction pp will be re-emitted in the flu¬ 
orescent channel and escape. Of the fraction pr of the 
photons re-emitted at resonance, a faction 1 — ft will be 
absorbed again before they are able to escape the shell. 
Of these [pr (1 — /?)], a fraction pp will be converted 
into fluorescent photons and escape [i.e., pfPr (1 — /3)]. 
Again, out of the resonantly re-emitted photons, a frac¬ 
tion 1 — /3 will be re-absorbed within the shell, contribute 
to the fluorescent re-emission and escape the outflow. It 
can be easily shown that, for each shell, the fraction of 
absorbed photons converted into fluorescent photons is 
given by: 


F f (t) = p F (1 - /3)] n , 


(18) 

v n—0 

red b 

while the fraction of absorbed photons that are able to 
escape will be: 


Fr(t) =pr/3^2[pr( 1- /3)} n . (19) 

n —0 

Because pr(1 — j3) < 1, the summation of the geometric 
series in Equations [T51 and fl9l converges, and 

F f =pf/[1-Pr(1-P)\ (20a) 

F R = pp R /[l-p R (l-p)\. (20b) 

Figure [5] shows the fraction of absorbed photons that 
escape in the fluorescent channel as a function of the shell 
optical depth for three representative values of p F . Anal¬ 
ogously to Eqns [13] and Qj] the blue and red resonant- 
emission components become: 


1 _ p -t(v) 

Fr(u) -x--d y 

max(x,l ) "iJ 


L 


,MSW = f 

J17 

ry°° 1 _ p~ T (y) 

(*) = / F R (y) --- 

Jvi 2 v 


em.red.MSV^yJ — 


d y, 


while the blue and red fluorescent components 
written as: 


(21a) 

(21b) 
can be 


-^em,blue,MS,F(^') 


*em,red,MS,F 


(x) 


fy°° 1 — p-rfy) 

/ F F (y) ---d y (22a) 

Jmax(x, 1) 

rv°° I _ p -t(v) 

/ F F (y) ---d y. (22b) 

Jyi 2 V 


5 In what follows, fluorescent transitions will be indicated with 
an *. 


In the right panel of Figure [4] we show the effect of ac¬ 
counting for multiple scatterings on the line profiles of the 
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-500 0 500 1000 1500 2000 -500 0 500 1000 1500 2000 

v„ bs [km/s] v obs [km/s] 

Fig. 4. — Left: Effect of fluorescent channel on the resonant P-Cygni profile. The Si” 1 " doublet line profiles are shown with and without 
the inclusion of the fluorescent emission (dashed and solid line respectively) for the single-scattering approximation. The pure absorption 
component of the profile is also shown for reference (dotted line). Right: Si"*” doublet profiles computed with single scattering approximation 
(black) and multiple scatterings (red). 




-400 -200 0 200 

V 0b . [km s"'] 


400 


Fig. 5. — The fraction of absorbed resonant photons re-emitted in 
the fluorescent channel after multiple scatterings depends on the 
gas column density as well as on the transition probability. We 
show the calculations for three transitions of Si”*”, as indicated in 
the label. 

Sill doublet. As expected, the re-scattered emission com¬ 
ponent at resonance is reduced significantly compared to 
the single-scattering approximation resulting in an in¬ 
creased intensity of the fluorescent lines. We note that 
this effect enhances the contamination of the Sill 1193A 
absorption component from re-emitted fluorescent pho¬ 
tons from the Sill 1190A transmission. This enhance¬ 
ment is particularly prevalent in low resolution spectra. 

2 .2. Spherical envelope observed with a circular finite 

aperture 

If the spectroscopic observations are made using an 
aperture that does not include the full extent of the scat¬ 
tering envelope, then the observed line profile can change 
dramatically. To illustrate the consequences, we consider 
here the case of a spherical envelope observed with a cir¬ 
cular aperture larger than the central source but smaller 
than the entire envelope, i.e., with R apel > Rsf and 
R aper < R\v ■ We also restrict our analysis to 7 = 1 . 

Clearly, the blueshifted absorption component of the 


Fig. 6 .— The P-Cygni line profile changes as function of size 
of spectroscopic aperture, from a pure absorption profile (/taper = 
.Rsf), to a classic P-Cyg profile (R ape r = Rw)- The profiles where 
computed assuming a spherical expanding envelope, 7 = 1 , r = 60, 
vo = 25 km s 1 and Voo = 450 km s _ 1 (we consider here a line 
with no fluorescent transition, such as, e.g., Silll A = 1260). 


profile will remain unchanged due to the presence of the 
aperture. The blue and red scattered emission, however, 
will change. In fact, for a shell of intrinsic velocity v, the 
aperture will block those photons scattered at velocities 
smaller than i> aper = v cos(9 ape r ), where 0 aper is such that 
sin[6 ape r ) = R ap er /r v . Clearly, see Figure [TJ each shell 
will correspond to a different 0 aper - 
In Figure [T| we show the edge-on view of a plane of 
constant v Q b s ( dotted line). As we saw earlier, only lay¬ 
ers with v > Vobs will contribute to the emission at ry>bs- 
Figure [T] shows that the effect of adding a circular aper¬ 
ture is to remove the contribution at u 0 b s from all shells 
with v > v up: where: 


Vuo = 


^obs 

COS0p>r)i 


(23) 


In terms of vq, and after a little algebra, we get: 
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2 2 
Vup = * 



(24) 


Thus, Iembiue^ x ) ( see Eq. HU) can be written here as: 


pv up 1 — p -t{v) 

IemMue(x) aper = / --- d y, (25) 

Jmax(x, 1) 

Analogously, I^ red {x) wiU be: 

rv up i _ P ~ T (v) 

We A{x) aper = - = -dy; (26) 

Jy i 2 2 / 

In Figure [ 6 ] we show how the P-Cygni profile changes 
with the ratio R apel /R sf- In the extreme case of 
Raper/R sf = 1 (dotted line), i.e., when the aperture 
is only as large as the source of continuum, the line is 
observed only in absorption, and blueshifted relative to 
the systemic velocity of the galaxy. A component of 
scattered re-emitted photons coming from the absorb¬ 
ing material is contributing to the blue side of the line, 
but no re-emitted photons are detected on the red side, 
because of both the effect of the aperture, and the -often 
neglected- occultation by the galaxy. 

As the ratio -Raper/ Rsf increases, fewer photons scat¬ 
tered toward the observer are blocked by the aperture. 
As a result, the emission component of the profile (cen¬ 
tered at the systemic velocity) becomes more and more 
pronounced. The shape of the absorption profile also 
changes because of the increasing contribution of photons 
scattered by material moving toward the observer. It is 
also evident from Figure [ 6 ] that the velocity at maximum 
absorption shifts toward higher blueshifted velocities as 
the the contribution from scattered emission increases 
(i.e., as -R a per/ Rsf increases). 

3. APPLICATION TO REAL SPECTRA 

As an example of its flexibility, we apply the SALT 
model to resonant line profiles observed in a stacked 
spectrum of Lya emitting galaxies. We first summa¬ 
rize the data and the measurements (Section 13.11) and 
then discuss the properties of that outflow that can be 
inferred from the absorption line analysis performed with 
the SALT model. 

3.1. Data and measurements 

The average spectrum modeled in this section was 
created by stackin g Cosmic Origin Spectrograph (COS 
lOreen et alJl2012l! medium-resolution spectra of a sam¬ 
ple of 25 known z ~ 0.3 Lya emitters (|Deharveng et all 
120081 : ICowie et all 1201 Oil . The details of the data reduc¬ 
tion and spectral extraction are presented in Scarlata et 
al. (2014). For each galaxy, we have accurate redshift 
meas urements obtained from the Ha emission line pro¬ 
files (ICowie et alJl2011 ). 

To create the stacked spectrum, we first blueshift the 
observed spectra of individual galaxies into the rest- 
frame using the measured Ha velocities. Then, at each 
wavelength we compute a flux-weighted average and a 
standard deviation. The mean stacked UV spectrum 
(shown after a box-car smoothing of 0.85A) is shown 
in Figure [71 where the shaded gray area corresponds to 


± one weighted standard deviation. In Figure [7] the 
top panel shows the number of galaxies that entered the 
stack at each wavelength. In the stacked spectrum we are 
able to clearly identify and reliably measure the features 
presented in Tables [2] and [3] and marked in Figure [7] 
The list includes five absorption lines and four fluores¬ 
cent emission lines. The list also includes the C ill A1175 
absorption line, which is mainly produced in the stellar 
photospheres. 

If we assume that lines are pure absorption and pure 
emission we could measure the bulk velocity of the gas 
from the peak velocity of the lines. We derive the peak 
positions by fitting Gaussian line profiles to the observed 
absorption/emission lines. When two lines of a given 
multiplet/ion are blended, we fit them simultaneously 
constraining the width of the Gaussian function to be 
the same for both lines. In Figure [5] we zoom-in on the 
spectral regions around different transitions in the Si + 
and Si ++ ions and plot the resulting best-fit Gaussian 
models. The errors on the peak wavelength were com¬ 
puted with a Monte Carlo simulation. We created 1000 
realizations of the stacked spectrum by changing the flux 
at each wavelength within ±1 a. The new profiles were 
fitted with a Gaussian, and the error on the peak wave¬ 
length was computed as the standard deviation of the 
1000 best-fit peak wavelengths. 

In Tables [2] and [3] we report the vacuum wavelength 
of the considered transitions, the observed peak wave¬ 
length of the profiles, as well as the velocity shift be¬ 
tween the galaxy’s rest frame velocity (computed from 
the Ha) and the peak velocity of the profiles. The stel¬ 
lar C III velocity is consistent with the systemic veloc¬ 
ity computed from the Ha emission line profiles. Note 
that the peak/trough velocities obtained from the Gaus¬ 
sian fits offer an easy mathematical representation of 
the data but do not add immediate physical meaning. 
The velocity profiles shown in Figure [ 8 ] are typical of 
star forming galaxies at both high and low redshifts (e.g . 

Shaplev et al Jl2003t iSteidel et al.l[20lltlJones et al.ll2012t 

Heckman et all 120111 : IWofford et al.l I2013T) . where they 
are usually interpreted as originating in a gas outflow 
probably driven by the current episode of star-formation. 

The results of the Gaussian fits presented in Fig¬ 
ure [5] would indicate that the gas is moving toward 
the observers with velocities -as measured at the max¬ 
imum absorption/emission- ranging between —160 and 
—220km s^ 1 . The average velocity computed in this way 
from all absorption lines is —185 ± 25 km s -1 . The pro¬ 
files in Figure [2] also show absorption at velocities as high 
as —500km s^ 1 , indicating the presence of multiple ve¬ 
locity components and/or a velocity gradient in the out¬ 
flowing gas. All detected fluorescent emission lines are 
blueshifted with respect to the systemic Ha velocity with 
peak velocities ranging between —88 and —137km s -1 . 
With an average outflow velocity of —100 ± 22km s -1 , 
the fluorescent emission components appear to have a 
systematically-lower velocity shift than the resonant ab¬ 
sorption lines. 

Resonant photons can be either re-emitted at reso¬ 
nance, or in the fluorescent channel, when available. 
Thus, the two lines originate in the same gas and will 
share the same kinematical properties. Naively, the mea¬ 
sured systematic difference in the bulk velocities of the 



















Scarlata & Panagia 



1200 1250 1300 1350 1400 

Wavelength [Angstroms] 


Fig. 7.— Composite rest-frame UV spectrum of 25 z ~ 0.3 Lya-emitting galaxies. Multiple absorption features are identified with vertical 
lines. Solid lines indicate stellar photospheric absorptions, dotted and dashed lines indicates resonant absorption originating in the galaxy’s 
ISM, with dotted and dashed showing low- and high- ionization metal lines, respectively. The top panel shows the number of galaxies 
used to compute the average spectrum at each wavelength. 
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absorption and fluorescent emission lines could then be 
interpreted as an indication that these lines formed in 
two kinematically-distinct components. In the following 
section, we use SALT to consistently model the scatter¬ 
ing from the outflowing gas, and show how the systematic 
velocity difference can be the result of a non-symmetric 
outflow. 

3.2. Modeling the line profiles 

Here we use SALT to model the observed absorp¬ 
tion/emission profiles observed in the stacked spectrum. 
The free parameters for the models are To, vq, and Uqo 
for the spherical outflow, and To, vo, foo and R a per , for 
the spherical outflow plus aperture. These parameters 
fully describe the density and velocity field of the galac¬ 
tic outflow, and do not depend on the particular transi¬ 
tion within a given ionic species. We therefore constrain 
the model’s free parameters by simultaneously fitting 
the four radiative transitions of Si + , observed around 
A = 1190A. We chose this spectral region because it 
shows the highest S/N of the stacked spectrum, and be¬ 
cause of the presence of both two resonant absorption 
and the corresponding fluorescent emission lines. To 
model the observed line profiles we use equation I22bl 
which accounts multiple scattering within each shell as 
explained in section [2] 

We derive the best-fit parameters for the symmetric 
outflow model with and without a view-limiting aper¬ 
ture, by performing a \ 2 minimization on the entire dou¬ 
blet profile, including both absorption and emission lines. 
The models were computed on the same velocity vector 
as the data, and then box-car smoothed with the same 
kernel, before proceeding to the computation and mini¬ 
mization of the x 2 . The best-fit line profiles for the Sill 
doublet are shown in Figure 0 and the parameter values 
are given in Table |4] The blue and red curves show the 
best fit for the spherical outflow with and without the 
spectroscopic aperture, respectively. 

The simplest spherical models (with or without aper¬ 
ture) well reproduce the depth, shape and central ve¬ 
locities of the blueshifted absorption components of the 
resonant doublet. This indicates that our simple assump¬ 
tions for the density and velocity fields of the scattering 
material are adequate representation of the Si + distribu¬ 
tion. On the other hand, the spherical model with no 
limiting aperture fails to reproduce the observed profiles 
in two key aspects: 1) it substantially overproduces the 
amount of scattered emission, and 2), due to the sym¬ 
metry in the considered configuration, it predicts that 
the emission component of the P-Cygni profile should 
be centered at the systemic velocity, while the observa¬ 
tions shows that the peaks of the fluorescent emission are 
clearly blueshifted. 

As we discussed in Section [2721 the effect of a spectro¬ 
scopic aperture is to selectively decrease the number of 
scattered photons that are able to reach the observer. As 
Figure [9] shows (blue curve), adding an aperture allevi¬ 
ates the first of the two discrepancies. However, because 
of the intrinsic symmetry of the outflow model, the scat¬ 
tered re-emission is still centered at zero velocity, and 
therefore the model still fails in fully reproducing the 
observed features. A blueshift in the scattered emission 
component originating in outflowing material can be ob¬ 
tained if the outflow is not spherically symmetric with 


respect to the central source. A simple way to imple¬ 
ment this, is by differentially weighting the contribution 
to the final profile from different portions of the enve¬ 
lope. Thus, we introduce a scaling factor - / 0 bsc - to 
the red component of the scattered profile (Ean fldl) : i.e., 
the radiation scattered in the half sphere moving away 
from the observer). Physically, the parameter f Q b sc can 
be used to mimic a face-on galaxy, where the disk is ab¬ 
sorbing part of the radiation emitted by the outflowing 
material (see Section |4j. The profile can now be written 
as: 

/(x) ^ = 1 Jabs T Ami,blue T /obsc X Ami,red- (27) 

/obsc, represents the fraction of Ami,red that is allowed to 
reach the observer. We refer to the above profile as the 
“asymmetric model”, to indicate that the receding half 
of the expanding envelope is seen less easily than the ap¬ 
proaching front part. The best fit asymmetric outflow 
model is shown in Figure [21 with the yellow line. This 
model is clearly able to simultaneously reproduce the rel¬ 
ative intensity of the fluorescent emission and resonant 
absorption lines, as well as their systematically different 
peak positions. 

We can test our results using a different resonant tran¬ 
sition in Si+. The resonant line at A = 1260.42 is ideal 
for this purpose. As Figure [8] shows, we detect both the 
resonant absorption and the corresponding fluorescent 
emission. This absorption originates in the same mate¬ 
rial where the Si IIA1190 doublet is produced and is thus 
perfect to test the parameters of the outflow model (i.e., 
the density, velocity field, and geometry). In Figure [TOl 
we show the observed Si IIA1260 profile, together with the 
model profile computed using the best-fit outflow param¬ 
eters derived from our analysis of the 1190-1193A dou¬ 
blet (i.e., changing only the transition dependent param¬ 
eters in the profile equations). Figure flOl shows that the 
model optimized to fit the Si II doublet fully reproduces 
the observed Si IIA1260 profile as well. In particular, it 
reproduces both the blueshifted absorption component 
as well as the intensity and peak wavelength of the fluo¬ 
rescent emission. We stress again that the parameters of 
the outflow are kept fixed to the best-fit values given in 
Table 01 

4. DISCUSSION 

Near and far-UV spectra include numerous resonant 
metal absorption lines that, combined with the appro¬ 
priate theoretical tools, provide powerful diagnostics for 
galactic outflows. High-quality rest-frame UV spec¬ 
tra are currently available for nearby individual galax¬ 
ies (e.g., with data from the Hubble Space Telescope), 
and stacked spectra of high-redshift galaxies (with data 
from 8m class telescopes). Soon, with the planned 30m 
telescopes, we will be able to study at high resolution 
absorption line spectra in individual objects up to the 
highest redshifts. 

Resonant blue-shifted UV absorption features are 
commonly modeled as originating in a thin shell 
of gas moving at the outflow velocity (as mea- 
sure d from the centroids of resonant absorption lines , 
e.g..lVerhamme et~a l. 2006; Sch aerer fc Verhammell200l : 

IVerhamme et al.l 120081 iSchaerer et al.l 1201 iff Substan¬ 

tial evidence, however, suggests that this description is 
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Fig. 8. — Zoom in on the spectral regions around the absorption and emission line features considered in this work. In each panel, the 
dashed vertical lines show the vacuum wavelength of each transition (as indicated by each label). The orange line shows the sum of the 
best-fit Gaussian profiles, while the dot-dashed line shows the continuum level. 



Fig. 9.— The absorption and re-emission profiles of the Si II dou¬ 
blet are well reproduced with an asymmetric model and accounting 
for the effect of the finite COS aperture size. Best fit model profiles 
to the Si II 1190-1193A doublet are shown for different outflow ge¬ 
ometries, as indicated in the label. All models shown in this Figure 
include multiple scattering within each shell. 


too si mplistic, as noted already by, e.g., iPettini et al.l 
(|2002f ). First of all, when the absorption lines are ob¬ 
served at high enough spectral resolution, they show 
asymmetric profiles covering a bro ad range of velocities 
(up to as much as -1 00 0 km s -1 . iTrcmonti et all 120071 : 
iDiamond-Stanic et al.ll2012t iSell et al.112014T) . This indi¬ 
cates that the gas is not confined to a thin shell, but 
rather is distributed in an extended envelope, with ve¬ 
locity and density changing with the distance from the 
galaxy. More realistic models of extended outflows, with 
velocity and density g r adients have been proposed (e.g., 
Prochakka et al.l 12011 : iRubin et al.1 1201 It IMartin et al.1 
201 .‘ill . and highlight the importance of properly account¬ 
ing for the geometry of the outflowing material. 

Second, well defined P-Cigny profiles from resonant 
transitions of Mgll, as well as the detection of fluo¬ 
rescent emission associated with resonant transitions of 
Si II, and Fell, are commonly observed dShaplev et al.1 
120031 : iFrance et al.1120101 IRubin et al.1120111) . indicating 


that the re-emission component from the outflowing gas 
cannot be neglected, particularly in compact galaxies or 
galaxies at high-z, where the spectroscopic aperture may 
include a substantial fraction of the extended scatter¬ 
ing outflow. Neglecting possible contribution from re¬ 
emitted photons may have important consequences for 
the determination of the gas column densit y and/or cov¬ 
ering f raction, as recently noted by, e.g..lProchask a et al.1 

(l20ll . 

The SALT model discussed in this paper provides a 
simple analytical description of the line profiles origi¬ 
nating in the expanding envelopes around galaxies, that 
properly account for multiple scattering of resonant pho¬ 
tons, scattered re-emission, and observational aperture 
effects. 

4.1. On the use of the absorption profiles as indication 
of covering fraction 

In the approximation of pure absorption and with well 
resolved line profiles, the gas apparent optical depth at 
a given velocity t(v) is often used to derive the appar¬ 
ent column d ensity profil e (—ln (I(v)/Io(v)) = t(v ) oc 
f\N(v), e.g. IPettini et all 1200211 . However, it is well 
known that the apparent column density obtained from 
line profiles can be underestimated if undetected satu¬ 
rated components contribute at some velocities. When 
two or more transitions of a given species differing only in 
the product of /A are available, then information about 
line saturation can be inferred from the comparison of 
the apparent optical depth profiles (that should be iden¬ 
tical within the observational errors). If saturation is 
present, the line with the highest oscillator strength will 
result in a lower apparent column density. The previous 
reasoning is correct only if the absorbing gas fully cov¬ 
ers the continuum source. If this is not the case, the line 
profile will also depend on the gas covering fraction (/c), 
as well as the optical depth (i.e., I/Io = 1 —/c(l —e _T )). 
Various works have used absorption-line profiles to deter¬ 
mine the gas fc , under the assumption that absorption 
lines are saturated (i.e., e~ T - > 0 in I/I 0 , e.g. I Jones et al.1 
120131 : IMartin fc Ik>uchdl 200 il) . When the scattering en¬ 
velope is included in the spectroscopic aperture, however, 
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this approach cannot be used: the resonant scattered 
re-emission affects substantially the residual intensity at 
line center, in different ways depending on the atomic 
level structure and the spontaneous transition probabil¬ 
ities.. 

We show this point in Figure 1111 where we present 
the resonant absorption profiles originating in an out¬ 
flowing envelope with f c = 1, for three transitions of 
Si + . The top panel shows the absorption component 
only of the profile: when only absorption is considered, 
the three profiles scales as expected, according to the 
value of /A for each transition. This is not true anymore 
when the resonant scattered and fluorescent re-emission 
are included (middle panel): close to the maximum ab¬ 
sorption, in fact, the transmission with the lowest value 
of /A is in fact the deepest! This is simply due to the rel¬ 
ative value of the branching ratios for the resonant and 
fluorescent channels for these particular transitions. At 
the largest velocities, where the filling from scattered ra¬ 
diation is less important and the optical depth is smaller, 
the absorption profiles are again proportional to /A. 

In the bottom panel of Figure [TT] we simulate the line 
profiles as they would be observed with COS, assuming 
a spectral resolution of 30km s _1 , and a noise of 10%. 
Because of the resolution and S/N ratio, the three lines 
are identical close to the core, and barely distinguish¬ 
able at the largest velocities. We note, however, that we 
did not account for uncertainties in the normalization of 
the spectra, which may affect the profiles, particularly 
at the largest velocities. The simulated lines also show 
that from the resonant profiles alone it is hard to iden¬ 
tify the presence of scattered re-emission: because at 
least 50% of the photons are re-emitted in the fluores¬ 
cent line, for all the transitions considered. However, the 
presence of the fluorescent emission is a clear sign-post 
that the absorption profiles will be affected by emission 
filling. Because the strength of the filling depends on the 
branching ratio for the fluorescent transitions, it is not 
correct to average together profiles of different lines of 
the same ion. Moreover, the effect of the filling of the 
absorption lines will also depend on the geometry of the 
outflow (see Section [3.21) as well as the the relative size of 
the scattering outflow region and the spectroscopic aper¬ 
ture.. 

As an example of what can be achieved by modeling 
absorption profiles with SALT, we have modeled the UV 
stacked spectrum of 25 Lya emitters in the nearby Uni¬ 
verse. The spectrum shows a number of resonant ab¬ 
sorption and fluorescent emission lines generating in the 
Si+ in the galaxy’s gaseous medium. The profiles of the 
resonant absorption lines are systematically blueshifted 
with respect to the systemic velocity, with an average 
gas velocity at maximum absorption of —185 ± 25 km 
s -1 . The peak velocity of the fluorescent emission lines, 
however, is systematically lower than the outflow veloc¬ 
ity obtained from the absorption profiles (—100 ± 22 km 
s -1 ). We showed with SALT how simple symmetric mod¬ 
els fail to reproduce this velocity difference, and asym¬ 
metry in the gas distribution needs to be present. Our 
best fit model with fobsc = 0.1 describes a geometry in 
which the half a sphere receding from the observer con¬ 
tributes only 10% to the fluorescent emission, which is 
then dominated by the emission from the half a sphere 


approaching the observer. Physically, this simple model 
can be used to describe the outflow from a disk galaxy 
seen approximately face-on, where the radiation scat¬ 
tered from the half-sphere receding from the observer 
has to go through the opaque disk, that will have an 
optical depth r = —/n(l — fobsc) = 2.3. 

The sample of galaxies that entered the stacked analy¬ 
sis is not randomly selected among starforming objects, 
but rather it comprises only galaxies showing Lya in 
emission. The conditions that allow more Lya radiation 
to escape from a galaxy are still highly debated, with 
some authors suggestin g that age is a dominant factor 
(e.g.- lCowie et al.ll201lD. while others advocating fo r dust 
extinction ('e.g. lHaves et al1l2013t lAtek et a1.1l2011l) . Our 
result seems to indicate that, on average, Lya-emitters 
are able to escape more easily in the direction perpen¬ 
dicular to the galaxy disk, than along the disk, suggest¬ 
ing that the viewing angle is also an important factor 
in determining the Lya escape/visibility. This result is 
not unexpected, since the direction perpendicular to the 
galactic disk is also the one that offers the minimum col¬ 
umn dens ity of diffuse material, a nd, thus, of neutral gas 
(e.g., see iVerhamme et al.ll2012l for results of radiative 
transfer simulations of Lya photons in a realistic spiral 
galaxy). In support of our findings, a recent structural 
analysis of hundreds of Lya emitters at 2 ~ 2.2 showed 
that Lya emitters tend to have smaller ellipticity than 
galaxies at similar redshift with no Lya emission, rein¬ 
forcing the idea that Lya photons escape more easily in 
the direction of minimum optical depth , i.e., minimum 
hydrogen column density (iShibuva et alJl20T3l . 
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Fig. 10. — Si IIA1260 absorption line profile and correspond¬ 
ing fluorescent emission, together with the models computed using 
the best-fit outflow parameters derived from the Si II 1190-1193A 
doublet. 

5. CONCLUSIONS 

In this work we have presented and discussed a simple 
Semi Analytical Line Transfer model, SALT, to describe 
the expected absorption and re-emission line profile gen¬ 
erated in a spherical outflow surrounding a galaxy with a 
finite size ( Rsf ) t under the Sobolev approximation. We 
derive the analytical profiles computed for the velocity 
field with the velocity increasing with distance from the 
galaxy (v oc r 7 ), and the densitity radial profile obtained 
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Fig. 11. — Neglecting emission filling of the absorption profiles 
can cause erroneous conclusions regarding the gas covering fraction. 
Top: absorption-only components of three resonant transitions of 
Si+, as indicated in the legend. The model correspond to a spher¬ 
ical envelope with vo = 35km/s, Uoo = 450km/s, and 7 ^ ^ = 30. 

When only the absorption components are considered, the profiles 
scale according to the value of /A. Middle: full line profiles includ¬ 
ing the scattered re—emission component. Particularly at velocities 
close to zero, the profiles do not scale as /A anymore. Bottom: 
simulated spectrum with 10% errors, and resolution of 30km/s. 

assuming a constant mass outflow rate (n;(r) oc (ur 2 )^ 1 ). 
We include the effect of multiple scatterings properly ac¬ 
counting for the atomic structure of the scattering ions. 
We also discuss how the line profile changes due to the 
effect of a circular spectroscopic aperture that does not 
cover the full extent of the outflow (Section [3]), and the 
case in which part of the outflow is obscured. 

Our analysis reproduces the main features observed in 
the profiles generated with more complex radiative trans¬ 
fer code s applied to gaseous ou tflows with the same ge¬ 
ometry ()Prochaska et alJ 12011 ). Namely, we show how 
-for a spherically symmetric outflow- the scattered re¬ 
emission is centered at zero velocity and thus alters the 
shape of the pure absorption component generated in the 
material in front of the emitting galaxy. Outflow veloci¬ 
ties computed from the wavelength of the absorption line 


trough do not include this effect. Such an analysis results 
in an overestimate of the outflow velocity. The intensity 
of the emission component of the line profile depends not 
only on the spectroscopic aperture used for the observa¬ 
tions, but also on the atomic structure of the particular 
ion used in the analysis and the spontaneous transition 
probabilities. In the case of the Si + 1190-1193A doublet, 
a photon absorbed in the 1190A transition has a higher 
chance of being re-emitted in the fluorescent transition 
to the 2P°3/2 level, thus somewhat limiting the “filling” 
effect on the absorption component. 

We have considered the resonance absorption and flu¬ 
orescent emission profiles observed in the average UV 
spectrum of 25 2 ~ 0.3 Lya emitters. With a simplistic 
Gaussian profile fit one would find that the average ve¬ 
locities computed from the trough (—185±25 km s 1 ) and 
the peak wavelengths (—100 ± 22 km s 1 ) of the profiles 
systematically differ, with the absorption-derived veloc¬ 
ity being more negative than the emission-derived veloc¬ 
ity. Regardless of the size of the spectroscopic aperture, 
a symmetric outflow of an arbitrary shape produces a 
profile with the emission component centered at the sys¬ 
temic velocity of the galaxy, and so cannot explain the 
systematic difference between the emission and absorp¬ 
tion velocities. On the other hand, we used SALT to 
show that this shifts comes naturally if most of the ra¬ 
diation scattered by the receding half of the outflow is 
obscured from view. This model can be interpreted as a 
simple representation of a disk galaxy observed face-on, 
where the thick disk acts as a semi-transparent screen 
for the backscattered radiation. This result thus indi¬ 
cates that, on average, galaxies tend to show Lya in 
emission more frequently when observed face-on. This 
idea is supported both by rece nt observations of high 
redshift Lya-emitting galaxies dShibuva et alj|20l4lh as 
well as high resolution radia tive transfer s im ulation of 
Lya photons in disk galaxies (jVerhamme et al1l2012f) . 

To conclude, by simultaneously reproducing both the 
resonant absorption and the associated resonant and 
fluorescent emission, the line profiles computed with 
our SALT model are more far-reaching than a simple 
absorption line based analysis. This is especially true 
when the data show evidence (e.g., the presence of a P- 
Cygni profile and/or fluorescent emission lines) of scat¬ 
tered re-emission from the galaxy wind. The formal¬ 
ism developed here can be easily extended to other ge¬ 
ometries, to account for dumpiness of the outflowing gas 
(Card et al., in prep), and to different ions/transitions. 
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TABLE 1 

Atomic data for Si II and Si III ions. Data taken from the NIST Atomic Spectra Databas^J 


Ion 

Vac. Wavelength A u i 

A s- 1 

Jlu 

Ei 

— E u 
eV 

9i ~ 9u 

Lower level 
Conf.,Term, J 

Upper level 
Conf.,Term, J 

Si II 

1190.42 

6.53x10 s 

2.77xl0 _1 

0 . 0 - 

10.41520 

2-4 

3s 2 3p 2P° 1/2 

3s3p 2 2P3/2 


1193.28 

2.69X10 9 

5.75xl0 _1 

0 . 0 - 

10.39012 

2-2 

3s 2 3p2P° 1/2 

3s3p 2 2P 1/2 


1194.50 

3.45X10 9 

7.37xl0 _1 

0.035613 

- 10.41520 

4-4 

3s 2 3p2P° 3/2 

3s3p 2 2P3/2 


1197.39 

1.40x10 s 

1.50xl0 _1 

0.035613 

- 10.39012 

4-2 

3s 2 3p2P° 3/2 

3s3p 2 2P 1/2 


1260.42 

2.57xl0 9 

1.22 

0.0 - 

9.836720 

2-4 

3s 2 3p2P° 1/2 

3s 2 3d2D 3/2 


1264.73 

3.04x10 s 

1.09 

0.035613 

- 9.838768 

4-6 

3s 2 3p2P° 3/2 

3s 2 3d2D 5/2 


a http://www.nist.gov/pml/data/asd.cfm 
b http://www.nist.gov/pml/data/asd.cfm 


TABLE 2 

Absorption lines measured in the stacked spectrum. 


Species 

Formation A va c 

A 

Aobs 

A 

An a 

Si ii 

ISM 

1190.42 

1189.83 ±0.09 

160 ± 25 

Si ii 

ISM 

1193.29 

1192.60 ±0.09 

174 ± 24 

Si ii 

ISM 

1260.42 

1259.56±0.08 

-218 ±21 

C ii 

ISM 

1334.53 

1333.81 ±0.10 

-164± 24 

Si hi 

ISM 

1206.50 

1205.75 ±0.08 

-188±20 

C hi 

Photo 

1175.53 

1175.42±0.1 

-28 ±25 


a Velocity shift of the absorption trough with respect to the 11 a 
emission line systemic velocity. 


TABLE 3 

Fluorescent emission lines measured in the 
stacked spectrum. 


Species 

Formation A V ac 

A 

Aobs 

A 

Av a 

Si ii* 

ISM 

1194.50 

1194.10 

- 100.2 

Si ii* 

ISM 

1197.39 

1197.04 

-88.2 

Si ii* 

ISM 

1265.00 

1264.42 

-137.5 

C ii* 

ISM 

1335.71 

1334.53 

-55.7 


a Velocity shift of the absorption trough with respect 
to the Hct emission line systemic velocity. 


TABLE 4 

Best fit values for the Si II lambdall90 — 1193 doublet. 


Model 

u 0 

[km s -1 ] 

Voo 

[km s -1 ] 

^0 

■^aper 

Rsf 

/obsc 

Spherical model full view 

38 

425 

160 



Spherical model limited view 

40 

426 

120 

2 


Asymmetric model - limited view 

55 

425 

45 

3 

0.1 
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